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ABSTRACT 

The radio polarization characteristics of milhsecond pulsars (MSPs) differ significantly 
from those of non-recycled pulsars. In particular, the position angle (PA) swings of 
many MSPs deviate from the S-shape predicted by the rotating vector model, even af- 
ter relativistic aberration is accounted for, indicating that they have non-dipolar mag- 
netic geometries, likely due to a history of accretion. Stokes tomography uses phase 
portraits of the Stokes parameters as a diagnostic tool to infer a pulsar's magnetic 
geometry and orientation. This paper applies Stokes tomography to MSPs, generaliz- 
ing the technique to handle interpulse emission. We present an atlas of look-up tables 
for the Stokes phase portraits and PA swings of MSPs with current-modified dipole 
fields, filled core and hollow cone beams, and two empirical linear polarization models. 
We compare our look-up tables to data from 15 MSPs and find that the Stokes phase 
portraits for a current-modified dipole approximately match several MSPs whose PA 
swings are flat or irregular and cannot be reconciled with the standard axisymmetric 
rotating vector model. PSR J1939-H2134 and PSR J0437-4715 are modelled in detail. 
The data from PSR J1939+2134 at 0.61 GHz can be fitted well with a current-modified 
dipole at (a, i) = (22±2°, 80±1°) and emission altitude 0.4 tlc- The fit is less accurate 
for PSR J1939-H2134 at 1.414GHz, and for PSR J0437-4715 at 1.44 GHz, indicating 
that these objects may have a more complicated magnetic field geometry, such as a 
localized surface anomaly or a polar magnetic mountain. 
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1 INTRODUCTION 

The two tools used most frequently to characterise the orien- 
tation and magnetic geometry of a radio pulsar are its pulse 
profile and position angle (PA) swing. T he rotating vec- 
tor model iRadhakris hnan fc Cookg|l969l ). which assumes 
an axisymmetric magnetic field, predicts an S-shaped swing 
across one pulse period and is traditionally used to deter- 
mine the inclinations of the magnetic axis of symmetry and 
the observer's line of sight to the rotation ajcis. However, hm- 
itations arise when analysing only the PA swing, especially 
as the magnetosphere is not axisymmetric in general, e.g. 
the magnetic field include s a current-modified component 
jHibschman fc Aro n^'200l''). 

In IChung fc'X lclatos ( 201(!) (hereafter CMIO), Stokes 
tomography was introduced as a diagnostic tool to be used 
alongside more traditional methods of analysis. It exploits 
the fact that the phase portraits traced out by the four 
Stokes parameters, when plotted against each other over 
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one pulse period, are unique for any given magnetic ge- 
ometry and orientation. An atlas of look-up tables, con- 
taining Stokes phase portraits and PA swings, was gener- 
ated by CMIO for a variety of simple models, including 
pure and current-modified dipole fields, filled core and hol- 
low cone beams, and the associated linear polarization pat- 
terns. CMIO also showed that, from a sample of 24 nominally 
"dipolar" pulsars, which obey the period-pulse-width rela- 
tion and/or exhibit clean S-shaped PA swings, the Stokes 
phase portraits of 16 objects are either inconsistent with 
low-altitude emission from a pure dipole field, or have a 
highly asymmetric surface emission pattern. 

In this paper, we turn our attention to millisecond 
pulsars (MSPs). Polarimetric studies of MSP radio emis- 
sion have uncovered complex behaviour not normally seen 
in slower pulsars. In particular, the PA swings of many 
MSPs arc neither c lean nor S -shapcd; instead, they are flat 
(Stairs ct al. 1999: lOrd et al.i 20 04). highly distorted (e.g. 
PSR J0437-4715; iNavarro et"ar.lll997l ) or ex tremely noisy 
(e.g. giant pulses from PSR J1824-2452A; iKnight et at] 
I2OO6I ). Additionally, the PA swing varies dramati cally with 
frequency for many MSPs l| Johnston et al.ll2008l ). indicat- 
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ing that the magnetic geometry changes a lot with altitude, 
or that the observed pulse profile comprises emission from 
several different regions and altitudes. 

The above trends suggest that MSPs have nondipo- 
lar magnetic fields. In a non-recycled pulsar, a dipole field 
can be distorted by several mechanisms, e.g. a current 
flowing alo ng the field lines (|Hibschman fc Arond I2001I : 
iDvkd 20081). or rotatio nal swccpback near the light cylin- 
der ( Hibschman fc Arons 2001 ; Dvks fc Harding 20ol ; iDvkd 
I2OOS ). In a recycled pulsar with a history of prolonged ac- 
cretion, another set of mechanisms comes into play. For ex- 
ample, accreted material channeled onto the magnetic poles 
distorts the fr ozen-in magnetic field as i t spreads towards 
the e quator llMelatos fc Phiiinevl I2001I: iPavne fc Melatod 



I2OO4I : Izhang fc Koiimal I2OO6I : IVigelius fc Melatod l20oir 
Quadrupolar magnetic fields, propo sed to expl ain the X- 
ray light curves of Her X-1 (Shakur a et aL 199l[). can eve n 
be comparable to the dipolar component I Long et al.ll200^ . 
Multipole fields can also be generated near the inner edge 
of t he partially diam agnetic accretion disk of an X-ray pul- 
sar l|Lai et al.ll 19991 ). Alternatively, as the pulsar is spun up 
by accretion, the magnetic pole drifts towards the rotation 
axis, dragged inwar d by the motion of superfluid vortices 
in the pulsar's core JSrinivasa n et al.l |l990l : lRudermanlll99ll : 
ICheng fc Dailll997l : lLamb et al. 2009). 

In this paper, we apply Stokes tomography to millisec- 
ond pulsar data drawn from the European Pulsar Network's 
(EPN) online database. In Section O we briefly review the 
fltting recipe for determining the optimal orientation and 
beam polarization patterns from observed pulse proflles and 
Stokes phase portraits. We also extend the model in CMIO to 
treat interpulse emission. We compare our improved look-up 
tables of Stokes phase portraits and PA swings to observa- 
tions of 15 MSPs in Section|3]to identify general trends. We 
then conduct detailed modelling of PSR J1939-I-2134, which 
has a strong interpulse, and PSR J0437— 4715, which has a 
pulse with multiple peaks, in Sections |4] and [5] respectively. 
Our flndings are summarised in Section [S] 



2 STOKES TOMOGRAPHY 
2.1 Radiation field 

For the convenience of the reader, we begin by summarising 
briefly how to determine the emission point and hence the 
polarization state of the radiation as a function of pulse 
longitude, following the recipe laid out in Section 2 of CMIO. 
Our notation and definitions copy CMIO. 

We define two reference frames, as in Figure 1 of CMIO: 
the inertial frame, in which the observer is at rest, with axes 
(ea;,ey,ez), and the body frame of the pulsar. The relative 
motion between the frames is computed by solving Euler's 
equations of motion (including precession in general but not 
in this paper). The line-of-sight vector w is chosen to lie in 
the By-Bz plane, making an angle i with e^. The rotation 
and magnetic axes lie along 62 and one of the body frame 
axes (es) respectively, separated by an angle a. We define 
a spherical polar grid (r, 6, (j>) in the body frame covering 
the region Xmin ^ f/^LC ^ O^f^TT, 0^0^ 

27r, with 64 x 256 x 128 grid cells, where the line 9 = 
lies along 63, and tlc = c/^ is the light cylinder radius. 



In this paper, we take Xmin = 0.01 and Xmax = 0.83 to 
accomodate the relatively small magnetospheres (and hence 
emission altitudes) of MSPs. 

Radiation from highly relativistic particles flowing 
along magnetospheric fleld lines is narrowly beamed. Hence, 
without relativistic aberration, the observed emission point 
xo(t) at any time t is located where the magnetic vector 
B[xQ(t),t] points along w. When aberration is included, 
the emission point xo( t ) at time t satisfies the equation 
(|Blaskiewicz et al.lll99ll : lDvkdl2008l l 



±t -I- f2 X xo/c 
|±t + n X xo/c 



(1) 



where t = B[xo(t), t]/|B[Xy(t), t] | is the unit tangent vector 
to the magnetic field at xo(t), and f2 is the angular veloc- 
ity vector. CMIO considered emission from only one pole for 
simplicity [i.e. -l-t in equation ([l])], and hence ignored in- 
terpulse emission. In this paper, we include emission from 
both the north and south poles, requiring both ±t terms 
to be retained in U]). At every instant, we thus have four 
emission points which satisfy ([T|), two in the hemisphere op- 
posite the observer (which he cannot see), and two facing 
the observer, which we label Pi and P2. We search the grid 
at a fixed altitude ro to find Pi and P2 ; the locations of Pi 
and P2 change with time in both the body frame and the in- 
ertial frame . The i7 x xp term in (flj) en codes the aberration 
effect, as in iHibschman fc Arond I2001I ). It is correct to or- 
der 0(r/r-Lc) and should be replaced by the full relativistic 
expression when warranted by confidence in the model and 
data. 

The Stokes parameters (/, Q, U, V) associated with the 
complex electric field vector E at xo(t), which describe the 
polarization state, are defined as 



/ = 



\Ey 



U 
V 



2Re{E^E*y) 
2lm{E^E;), 



(2) 
(3) 
(4) 
(5) 



where 7 is the polarised fraction of the total intensity, 
L = {Q^ + U'^y^^ is the linearly polarised component, and 
V is the circularly polarised component. The observed elec- 
tric field vector E, assumed to be in the direction of the 
particle acceleratiorQ, is the incoherent sum of the elec- 
tric field vectors at Pi and P2, viz. E = Ei -I- E2, where 
the relative phase between Ei and E2 fluctuates randomly. 
The observed Stokes parameters therefore reduce to 7 = 
7i -I- 72, Q = Qi -I- Q2, and U = Ui + 1/2. In this paper, we 
assume that all the emission is linearly polarised for simplic- 
ity, i.e. V = 0. Circular polarization will be examined in a 
companion paper. 

The X- and y- components are measured with respect 
to an orthonormal basis (x, y) which is flxed in the plane 
of the sky. In this paper, we choose x = f2p/|f2p| and y = 
X X w, where flp = ft — {CI ■ w)w is the projection of fl 
onto the sky. Then the polarization angle, 7/), between x and 



^ The instantaneous acceleration is inclined slightly with respect 
to the normal (or binormal) of B at xo(t), because the emit- 
ting charges corotate. For more details, see t he discussio n around 
equation (2) in CMIO and equation (A3) in lDvksl 1I2OO8I) . 
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the linearly polarised part of E is given by 
V. = -tan -. 



(6) 



The observational data obtained from the EPN are not 
necessarily expressed in the canonical basis (x, y) . However, 
the Q-U phase portrait has the same shape in any Carte- 
sian basis; if the basis is rotated by an angle /3 with respect 
to X and y, the Q-U phase portrait rotates by an angle 
2/3 without being distorted, unlike the I-Q and I-U phase 
portraits, which change shape. Hence, when analysing the 
data, the first step is to reproduce the shape of the Q-U 
phase portrait as closely as possible without worrying about 
the orientation; infer j3; rotate the (/, Q, U, V) data into the 
canonical basis provisionally defined through /3; and then ad- 
just a, i, and the beam and polarization patterns iteratively 
to reproduce the I-Q and I-U portraits. The recipe for doing 
so is explained in Section 2.6 and Figure 3 of CMIO. 



2.2 Look-up tables of Stokes phase portraits 



Figures I18H33I in the Appendix display look-up tables of 
Stokes phase portraits and PA swings, similar to those in 
CMIO, updated to include interpulse emission. The figures 
are organised into four groups, corresponding to two beam 
models (filled core and hollow cone) and two polarization 
models {L oc cos 8, L oc sm9\ see CMIO). All the look-up ta- 
bles are for a current-modified dipole magnetic field (CMIO) 
composed of a pure dipole plus a toroidal component with 
magnitude 

-B^ = — -BpCosasin6'r/r-LC, (7) 
where Bp = {B^ -\- Bg)^^^ is the poloidal field strength. 



2.3 Inter pulses 

An interpulse is a secondary pulse separated from the 
main pulse by approxim ately 180° of rotational phase 
iManchester fc Lvnelll977l ). It is believed to arise when a 
pulsar is a nearly orthogonal rotator viewed nearly side-on, 
shining from both magnetic poles, i.e. with a ~ i ^ 90°, 
where 'R i' means 'with in roughly one beam width' in this 
context (|Petroval |200i). 

Figure[T]compares the Stokes phase portraits, pulse pro- 
files, and PA swings for a pure dipole and current-modified 
dipole emitting from one and two poles for one illustrative 
orientation {a,i) = (80°, 70°). For clarity, relativistic aber- 
ration is not included in this example (compare Section [2.41 
et seq.). Clockwise from the top left panel, the figure displays 
(i) a pure dipole with no interpulse, (ii) a current-modified 
dipole at r = 0.13rLC with no interpulse, (iii) a pure dipole 
with an interpulse, and (iv) a current-modified dipole at 
r = 0.13rLC with an interpulse. The top two panels in Fig- 
ure [1] which have no interpulse, are the same as in Figures 
5-8 and 30-33 in the look-up tables in CMIO. 

Figure [2] shows the loci xo(i) traced out by Pi and P2 
over one rotation in the body frame of the pulsar for the 
various cases in Figure [1] The panels are arranged as in Fig- 
ure[T] The bottom panels, in which the interpulse is present, 
show two paths, one in the north hemisphere, and one in the 
south. For the current-modified dipole (right panels), the 




V 
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Figure 1. Stokes tomography of a model pulsar with an inter- 
pulse but without relativistic aberration. Top left: dipole field, 
no interpulse. Top right: current-modified dipole emitting at r = 
O.lSrLCi no interpulse. Bottom left: dipole field with interpulse. 
Bottom right: current-modified dipole emitting at r = O.lSrLC 
with interpulse. Within each quadrant of the figure, the five sub- 
panels display (clockwise from top left): ///max and PA (in ra- 
dians) as functions of pulse longitude, I-Q, Q-U, and I-U. The 
orientation is («,«) = (80°, 70°). The beam pattern is given by 



loci are asymmetric, as discussed in CMIO. For definiteness, 
we consider a fiUed-core beam, viz. 

-1/2 r 
jexp 



i{e. 



= (271 



V(2a^)] 



+exp[-(e-^)V(2<T^)]}: 



(8) 



which is represented by greyscale shading in Figure [2] In (|8]), 
a is the width of the beam, chosen arbitrarily to equal 10°. 
We also choose the linear polarization pattern to be 



L{e,(p) = i{e,(p)\cosi 



(9) 



in Figures [T] and [21 Other choices (e.g. L oc sinS) are equally 
valid and have been found empirically by CMIO to match 
the observational data in many objects. 

In Figure [TJ the interpulse traces out a small, secondary 
loop within the primary pattern in the I-Q, I-U and Q-U 
phase portraits. It also changes slightly the range of U and 
Q covered by the main pulse. For example, for the dipole 
(left panels), the maximum value of U decreases from 0.9 
to 0.8 with the addition of the interpulse. The size of the 
secondary loop (i.e. the intensity of the interpulse) increases 
relative to the main pulse as a approaches i, as expected. 
Along a = 90° and i = 90°, when the interpulse and the 
main pulse peak at the same intensity, the primary and sec- 
ondary patterns overlap in the I-Q, I-U and Q-U planes, 
and the phase portraits are indistinguishable from the non- 
interpulse case. The shapes do not overlap exactly for other 
orientations, where the main pulse is brighter than the in- 
terpulse. 

In the I-U and Q-U planes, the balloons and heart 
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Figure 2. Locus xo(t) traced out by the emission point(s) P\ and 
P2 over one rotation, in the body frame of the model pulsar con- 
sidered in Figure[T] Top left: dipole field, no interpulse. Top right: 
current-modified dipole emitting at r = 0.13rLCi no interpulse. 
Bottom left: dipole field with interpulse. Bottom right: current- 
modified dipole emitting at r = O.lSrLC with interpulse. The 
beam pattern 1(9, <f>) is represented by greyscalc shading (bright- 
ness oc /). 

shapes seen in CMIO are also seen when an interpulse is 
present. For example, for a = i = 90°, the Stokes param- 
eters trace out two reflection-symmetric patterns with pos- 
itive and negative U to form complex, interlocking shapes 
(see Figures I18H32I from the atlas of look-up tables in the 
Appendix) . As expected, the patterns are more intricate for 
a hollow cone than for a filled core. For example, the Q-U 
portrait at (a,i) — (70°, 80°) for a filled core contains an 
asymmetric, tilted heart shape and a small oval, both con- 
nected at Q = f/ = (Figure [20)) . The same orientation for 
a hollow cone shows a broader heart shape with two large, 
secondary ovals (Figure [2HJ- 

2.4 Relativistic aberration 

In the observer's reference frame, charged particles flowing 
outwards ultra-relativistically along poloidal magnetic fleld 
lines also have a small transverse velocity component be- 
cause they corotate with the star as part of the highly con- 
ducting magnetosphere. This displaces xo(t) by a distance 
of order r/rLC compared to its position when aberration 
is neglected. The electric field vector (parallel to the parti- 
cle's acceleration vector) is also displaced, re sulting in the 
well-known delay-radius relation jBla skiewicz et al.l Il99ll : 
iHibschman fc Arond [2OO1I : iDvksl l2008i ). According to this 
relation, the centre of the pulse profile leads the steepest 
point of the PA swing by 4r/rLC- 

We compute xo(t) directly, including aberration, by 
solving ([1]) numerically. As a cross check, we compare the 
numerical solution with the analytic approxim ation given 
by equation (F2) of IHibschman fc AronsI ()200lh . where the 
tangent field at the aberration-shifted emission point, t, can 
be expressed as the sum of the tangent field at the origi- 



nal, non-aberrated emission point, to, plus a perturbation 
ti = (f2 X xo)/c — [{ft X xo)/c ■ to]to. As the aberration- 
induced deflection angle grows linearly with r, equation 
(F2) holds most accurately for small r. At r = O.lSrLC 
and O.SlrLC, the direct and approximate calculations of t 
agree to within ^^10% and ~20% respectively. Note that, 
although we calculate t directly from fl]), the expression 
±t -I- f2 X xo/c in ([ij itself breaks down near the light cylin- 
der, where quadratic relativistic corrections come into play. 

Figure [3] illustrates how aberration modifies the Stokes 
phase portraits, pulse profiles, and PA swings for pure dipole 
and current-modified dipole magnetospheres. For the sake 
of clarity, we do not include interpulse emission in Figure 
[S] although, in general, interpulse and aberration effects are 
additive, as one can tell from the look-up tables in the Ap- 
pendix. The emission is placed arbitrarily at an altitude of 
O.lr-LC to ensure a reasonably strong effect. We note that 
aberration introduces an altitude dependence in the case of 
a pure dipole, which is absent in the non-aberrated dipole 
considered in CMIO. 

Aberration acts mainly to shift the relative phases of 
the pulse centroid and PA swing inflection point. To lowest 
order in r/rLC, the radius-delay relation predicts that the 
pulse profile is phase shifted by « —r/r-^c radians, whereas 
the PA swing is phase shifted by ~ 3r/rLC radians. Figure 
[3] shows that, for r — O.lrLC, the pulse profile is shifted 
by ~ —0.08 radians, whereas the PA swing is shifted by 
« -1-0.27 radians. These shifts are enough to dramatically 
broaden the I-Q pattern, twist the I-U pattern, and tilt the 
Q-U pattern for the pure dipole (left panels of Figure [S]). For 
the current- modified dipole (right panels), the I-Q pattern 
narrows, the I-U patterns twists and rotates, and the Q- 
U pattern rotates. Figure |3] shows the loci of xo(t) traced 
out by Pi over one rotation (there is only one set of emission 
points without an interpulse) , with each panel corresponding 
to the cases in Figure |31 The loci of the aberrated emission 
points (bottom panels) are shifted in (f) relative to the non- 
aberrated points (top panels). 

In CMIO, it is shown that, for a pure dipole field with- 
out aberration, all phase portraits are reflection symmetric 
about U — 0. Aberration breaks this symmetry, causing the 
shapes in the I-U and Q-U plane to tilt (see Section 4 in 
CMIO). Aberration also changes the tilt and the relative 
sizes of the shapes in the phase portraits, e.g. the ventricles 
of the hearts in the U -Q plane. 

2.5 Tilted axis of symmetry 

If the beam pattern is centred on another axis that is slightly 
tilted with respect to the magnetic axis (e.g . fan beams in 
the o uter magnetosphere; ICheng et al.ll2000l : IWatters et al.l 
|2009| ). the pulse profile is also phase-shifted relative to the 
PA swing. Tilting the beam axis away from the magnetic 
axis can therefore mimic closely (though not exactly) the 
effects of aberration (see Section 2.4.1 in CMIO). To illus- 
trate. Figure [5] compares the phase-shift caused by the tilt 
of the beam axis to that caused by aberration. We plot six 
pulse profiles and PA swings for a pure dipole with (q, i) = 
(30°, 40°). The three panels on the left correspond to beam 
axes which are tilted with respect to the magnetic axis by 
(6l',<^') = (10°, 0°) (top), (10°, 45°) (middle), and (10°, 90°) 
(bottom), all emitting at r = 0.02rLC- Using the top panel 
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Figure 3. Stokes tomography for a model pulsar including rel- 
ativistic aberration but without an interpulse. Top left: dipole 
field, no aberration. Top right: current-modified dipole emitting 
at r = O.frLC) no aberration. Bottom left: dipole field emitting 
at r = O.frLC with aberration. Bottom right: current-modified 
dipole emitting at r = O.frLC with aberration. Within each quad- 
rant of the figure, the five subpanels displays (clockwise from top 
left): ///max and PA swing (in radians) as a function of pulse 
longitude, I-Q, Q-U , I-U. The orientation is {a,i) = (30°, 40°). 
The beam pattern is given by JHJ. 
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Figure 4. Locus xo(i) traced out by the emission point Pi across 
one rotation, in the body frame of the model pulsar considered in 
Figure|3] Top left: dipole field, no aberration. Top right: current- 
modified dipole emitting at r = O.lrLCi aberration. Bottom 
left: dipole field emitting at r = O.lrLc with aberration. Bot- 
tom right: current-modified dipole emitting at r = O.frLC with 
aberration. The beam pattern /(6, 0) is represented by greyscale 
shading (brightness oc /) . 
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Figure 5. Pulse profiles and PA swings for pure dipoles with 
identical orientations (a, i) = (30°, 40°). The left panels show 
beam axes offset by (0',0') with respect to the magnetic axis, 
corresponding to (e',(^') = (fO°,0°) (top), (9',^') = (f0°,45°) 
(middle), and {9',<f>') = (10°, 90°) (bottom). In the left panels, 
the emission radius is fixed at r = 0.02rLC- The right panels cor- 
respond to emission heights r = 0.02rLC (top), O.OTrLC (middle), 
and O.frLC (bottom). In the right panels, the magnetic and beam 
axes are aligned. 



as a reference point, the pulse centroid leads the PA swing 
by — O.lSrad in the middle panel, and by — 0.36rad in the 
bottom panel. On the right-hand side of Figure [5] the beam 
axis and magnetic axis are aligned, but we vary the emission 
radius. As aberration causes the pulse centroid to lead the 
PA swing by — 4r/rLC, one could just as well attribute the 
phase shifts seen in the left-hand side to the emission radius 
increasing from r = 0.02rLc (top), 0.07rLC (middle), and 
O.lr-LC (bottom). 

As it is not possible to distinguish the two effects with- 
out additional a priori information, throughout this paper 
we assume the beam is centred on the magnetic axis, except 
in Sections [3] and [S] where we work with tilted beams re- 
constructed empirically from the pulse shape (because the 
magnetic-pole-centred model does not fit the data). 



3 MINIATURE POPULATION STUDY 

We now survey the pulse profiles, Stokes phase portraits, and 
PA swings for a selection of 16 MSPs from the EPN online 
databas^ (|Lorimer et al.|[l998l ). These objects are chosen 
because they have pulse periods < 10 ms, except for PSR 
J1022-I-1001, which was included because of its interesting 
Stokes phase portraits and because its orientation angles a 
and i have been measured with some degree of confidence 



^ Available at: http:/ /www.mpifr-bonn.mpg.de/pulsar/data/| 
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by previo us authors (IStairs et al ] | 19991 ). All objects were ob- 
served bvlStairs et all (Il999l) except for PSR J0437-4715 
(jManchester fc Johnstonll 19951 ). We note that the PA swing 
lb defi ned in this paper corresponds to —ip in IStairs et al.l 
l|l999h . Table [3] quotes the size of the magnetosphere for 
each MSP (in units of the stellar radius, r*) and the frequen- 
cies where EPN data are available. In the few cases where 
rotating vector model fits have been attempted in the liter- 
ature, a and i are also quoted, together with the published 
uncertainties. 



3.1 General trends 

In Figures [BHSl we present the Stokes phase portraits, pulse 
profiles, and PA swings for the eight pulsars with the clean- 
est data. All Stokes parameters are normalized by the peak 
intensity /max. MSPs generally have a lower degree of lin- 
ear polarization than non-recycled pulsars, so their phase 
portraits are correspondingly noisier. The PA swing is only 
drawn at pulse longitudes satisfying L ^ O.l-Lmax and 
/ ^ 0.1/max, where Lmax is the peak value of L. As the 
absolute orientation of flp (and hence the angle /3 between 
the measured and canonical bases) for each set of data is 
unknown, we start by analysing just the shape of the Q-U 
phase portraits, as discussed in Section [2. II and CMIO. 

Many of the pulse and linear polarization profiles are 
highly asymmetric, suggesting a complex emission pattern. 
In objects where I{t) and L[t) have multiple peaks or inter- 
pulse emission, the Stokes phase portraits feature multiple 
loops, each corresponding to an individual peak. For ex- 
ample, each of the five peaks seen in PSR J0437— 4715 at 
1.44 GHz (Figure [6l top row) corresponds to a distinct sub- 
pattern in the I-Q and I-U planes (see Section [5|. In the 
Q-U plane, the pattern formed is an asymmetric figure-eight 
over a slightly curved line. PSR J1022-I-1001 (Figure|6] rows 
4-6) has an asymmetric, double-peaked pulse profile, which 
produces an asymmetric heart shape in the Q-U plane. Also 
interesting is PSR J1939+2134 (Figure H rows 4-5), which 
has a strong interpulse, whose phase portraits narrow with 
increasing frequency, while those of the main pulse broaden. 

The PA swings for the MSPs featured in Figures [BHH] are 
less informative. In several cases, where the PA swing is fiat 
or noisy, the phase portraits still trace out a recognisable 
pattern. For example, in Figure [G] the PA swing of PSR 
J0437— 4715 is fiat overall but punctuated by several dips, 
while the PA swing of PSR J1012-f 5307 is completely fiat. 
Their phase portraits, in contrast, reveal balloons and figure- 
eights. In FigureEl the PA swings of PSR J1713-I-0747, PSR 
J1744-1134 and PSR J1911-1114 are fiat, yet their phase 
portraits are distinguished by straight lines and balloons. 
The same is true of PSR J1939-I-2134 in Figure H This is 
another instance, to be added to those in CMIO, where the 
Stokes phase portraits carry important extra information 
which is not apparent from the pulse profile and PA swing 
alone. 



3.2 Magnetic geometry and orientation 

The data in Figures |SHH] are too low in quality to allow de- 
tailed fits for the angles a and i and the magnetic geometry, 
except for PSR J1939-f2134 and PSR J0437-4715, which 



we model in detail in Sections |4] and [5] It is still instructive, 
however, to compare the observed Stokes phase portraits in 
Figures [SHS] with the atlas of look-up tables in the Appendix 
and make some general remarks. 

With the exception of PSR J1022-I-1001, which has an 
S-shaped PA swing, the PA swings of the five other MSPs 
are flat or noisy, ruling out a purely dipolar magnetic ge- 
ometry. Below, we list the MSPs where we have been able 
to find approximately matching orientations, beam patterns 
and linear polarization models for a current-modified dipole. 
The relevant look-up tables are appended in parentheses. In 
inferring the orientations, we follow some rules of thumb, (i) 
If there is an interpulse present, we limit the look-up range 
to a « i > 60°. (ii) If there is more than one peak, we 
model the emission as a hollow cone, (iii) If L peaks with /, 
we assume L = I cos 8, whereas, if L vanishes at the pulse 
centroid, we assume L = /sin 6'. Away from the a = i di- 
agonal, both linear polarization patterns yield similar phase 
portraits. 

(i) PSR J1012+5307 (0.61 GHz, FigureEl thud row): hol- 
low cone, L = /cos6l, {a,i) = (70°, 10°) (Figures [2M29)l . 
This object has an interpulse. The balloons in the I-Q, I-U, 
and Q-U planes match approximately the orientations of the 
balloons in the model, although they have different widths. 

(ii) PSR J1713-h0747 (1.414 GHz, Figure [T] third row): 
filled core beam, L = IsinO, {a,i) = (80°, 30°) (Figures 
ll8H2ip . The straight lines in all three phase portraits match 
the model, although the gradient of the I-Q line is less steep 
than in the model. 

(iu) J1744-1134 (0.61GHz, Figure fourth row): filled 
core beam, L = /cos6i, (a,i) = (30°, 20°) (Figures [T8l42T|) . 
The straight line in the I-Q plane matches a thin balloon 
in the model, while the balloons in the I-U and Q-U planes 
match balloons in the model. Note that the balloons in the 
model are tilted upwards {dU/dQ > 0), whereas in the data 
they are tilted downwards {dU/dQ < 0). 

(iv) J1824-2452 (0.61 GHz, Figure 13 fifth row): hollow 
cone beam, L = I cos 9, {a,i) = (60°, 60°) (Figures [2M29)) . 
The pulse profile has three peaks, suggesting that this object 
may have a double-peaked interpulse. The data matches the 
model if the Q-U pattern is rotated by « 180°. 

Where multi- frequency observations are available, we 
only analyse the frequency at which the phase portraits are 
resolved best. We emphasize that the matches are approx- 
imate, and that the figures in the Appendix show only the 
phase portraits at one altitude, viz. r — O.lrLC. More de- 
tailed modelling of / and L as a function of emission altitude 
and {8, (j)) must be done to obtain more accurate matches, 
including the po ssibility that the emission originates from 
several altitudes l|johnston et al. II2OO8). 

For PSR J1022-I-1001 (Figure H rows 4-6), whose S- 
shaped PA swing is nominally dipolar, the heart shape in 
the Q-U plane roughly matches a pure dipole at {a, i) « 
(70°, 20°) for a hollow cone with either polarization model. 
However, the observed heart shape differs slightly from 
the model, and L{t) is actually triple-peaked, not double- 
peaked. Further information on the polarization basis (e.g. 
at several frequencies) is required in order to accurately de- 
termine the orientation and magnetic geometry. 

We now test whether the published a and i values in 
Table |3l inferred from the rotating vector model, are con- 
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Table 1. Pulse periods P and observation frequencies of 15 millisecond pulsars from the EPN online database with P < 10 ms (except for 
J1022+1001, which is included because of its interesting Stokes phase portraits). Where a rotating vector model fit has been published 
previously, a and i are quoted with their published uncertainties. The uncertainties for a a nd i for PSR J0437— 4715 are not given in 
the literature. All pulsars were observed bv iStairs et al.i (1999,), except for PSR J0437— 4715 jManchester fc Johnstoiill 19951) . 
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sistent with the observed Stokes phase portraits. We refer 
the reader to the look-up tables in Figures [1811321 in the Ap- 
pendix. For PSR J1022+1001 (Figure [1 rows 5-6), the PA 
swings at 0.61 GHz and 1.414 GHz imply two very different 
orientation s, namely, (a, i) = (140°, 135°) and (83°, 76°) re- 
spectively ijStairs et al.lll999l ). Already, this is worrying, as 
the orientation of a given pulsar should be unique, no mat- 
ter what altitude the emission comes from. Moreover, nei- 
ther of these orientations yield Stokes phase portraits which 
match the data, for any beam or linear polarization pattern. 
One can verify this easily by examining the phase portraits 
in the vicinity of (a,j) = (40°, 50°) and (80°, 70°) in Fig- 
ures [TSH321 in the Appendix. For example, for a hollow cone 
with L = I sin 6, at (80°, 70°), there are three interlock- 
ing ovals in Q-U, unlike the heart shape in Figure [S] For 
PSR J1824-2452 (Figure [3 fifth row), the rotating vector 
model predicts (a, i) = (41°, 81°). For a hollow cone, the 
look-up tables in Figure I26H28I show interlocking ovals in 
the Q-U plane, whereas the data reveal an oval joined to 
a straight line. Significantly, all these discrepancies are in 
the shape, not the orientation of the Q-U portrait, which is 
basis- independent (CMIO). We comment on the published 
{a,i) fits for J0437— 4715 in Section O where we model the 
object in detail. 



3.3 Emission altitude 

We now discuss how the pulse profiles and Stokes phase 
portraits evolve with frequency for PSR J1022+1001 and 



PSR J1939+2134, the only EPN MSPs with adequate muhi- 
frequency data. 

At 0.41GHz, PSR J1022+1001 has a double-peaked 
intensity profile. The first intensity peak is itself double- 
peaked in terms of its linear polarization, resulting in three 
L peaks overall. In the phase portraits, the first intensity 
peak corresponds to the bottom loop of the figure-eight in 
the I-Q plane and the large balloon in the I-U plane. At 
0.61 GHz, the second pulse is stronger than the first. In the 
I-U plane, the second pulse corresponds to the long, straight 
tail emerging from the bottom of the balloon. 

At 0.61 GHz, PSR J1939+2134 displays a single-peaked 
interpulse which peaks at « 0.6/max- The main pulse is also 
single peaked. The two pulses trace out qualitatively similar 
patterns on the I-Q and I-U planes, namely elongated bal- 
loons, whose major axes are tilted by ~ 0° and f» 20° relative 
to the Q = and U — axes respectively. At 1.414 GHz, 
the peak of the interpulse drops to ~ 0.4Jinax, and the main 
pulse is double-peaked. There is a dramatic difference in Q 
for the main pulse: the slope of the major axis of the bal- 
loon changes sign, from dQ/dl < (0.61 GHz) to dQ/dl > 
(1.414 GHz). The balloons of the main pulse also broaden, 
while those of the interpulse narrow. The second peak of the 
main pulse appears in the phase portraits as a kink in the 
I-Q and I-U balloons. 

For all the MSPs, the different ways in which individ- 
ual peaks evolve with frequency imply that I and L depend 
on 9 and in a complicated way. The profile components 
might originate from different emission regions whose mag- 
netic geometries are different functions of r. As the aberra- 
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Figure 6. Pulse profiles, Stokes phase portraits and PA swings for PSR J0437— 4715 at 1.44 GHz and 4.6 GHz llManchester fc JohnstonI 
1 19951) . PSR J1012+5307 at 0.61 GHz, and PSR J1022+1001 at 0.41 GHz, 0.61 GHz, and 1.414 GHz jStairs et al.lll999h . The data for 
each pulsar occupy a row in landscape layout. Prom left to right, the columns show: (1) ///max (solid curve) and ////max (dashed curve) 
versus time (in s), (2) the I-Q phase portrait, (3) the I-U phase portrait, (4) the Q-U phase portrait, (5) the PA swing versus time (in 
s) (data points with L ^ O.lLmax and / ^ 0.1/max plotted only). Data are presented courtesy of the EPN online archive. 
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Figure 7. Pulse profiles, Stokes phase portraits and PA swings for PSR J1713+0747 a t 0.41 GHz 0.61 G Hz and 1.414 GHz, PSR 
J1744-1134 at 0.61 GHz, PSR J1824-2452 at 0.61 GHz, and PSR J1911-1114 at 0.41 GHz dStairs et al.lll999h . The data for each pulsar 
occupy a row in landscape mode. From left to right, the columns show: (1) ///max (solid curve) and ////max (dashed curve) versus time 
(in s), (2) the I-Q phase portrait, (3) the I-U phase portrait, (4) the Q-U phase portrait, (5) the PA swing versus time (in s) (data 
points with L ^ O.lLmax and / ^ 0.1/max plotted only). Data are presented courtesy of the EPN online archive. 
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Figure 8. Puls e profiles. Stokes phase portraits and PA swings for PSR J1911-1114 at 0.61 GHz, and PSR J1939+2134 at 0.61 GHz 
and 1.414 GHz llStairs et al.lll999l) . The data for each pulsar occupy a row in landscape mode. Prom left to right, the columns show: (1) 
I/hwz.^ (solid curve) and hjl-^^y^ (dashed curve) versus time (in s), (2) the I-Q phase portrait, (3) the 1-V phase portrait, (4) the Q-V 
phase portrait, (5) the PA swing versus time (in s) (data points with L > O.lLmax and I > 0.1/max plotted only). Data are presented 
courtesy of the EPN online archive. 
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tion and toroidal field increase with r, they also distort the 
path xo(i), further complicating I{t) and L(t). 



4 A DETAILED EXAMPLE OF INTERPULSE 
EMISSION: PSR J1939+2134 

In this section and the next, we model the pulse and 
hnear polarization profiles of PSR J1939+2134 and PSR 
J0437— 4715 in detail and attempt to determine their ori- 
entation from their Stokes phase portraits. 

We first apply the iterative recipe from CMIO to 
PSR J1939+2134. which has P = 1.558 ms and P = 
1.051 X 10"^®s s"^ (|Kaspi et al.ll 19941 ') . making it the second- 
fastest known MSP. Data for this object, at 0.61 GHz and 
1.414 GHz, are obtained from the EPN online archive . The 
data were originally published in IStairs et al.1 (|l999h . We 
choose this object because of its strong interpulse emission. 
As mentioned in CMIO, the data published in the EPN are 
not expressed in the canonical polarization basis described 
in Section [21 Additionally, the emission altitude for this ob- 
ject has not been estimated at either frequency due to the 
flatness of its PA swing. We are therefore obliged to infer /? 
indirectly, from qualitative considerations, in order to bring 
the data of both frequencies into the canonical basis. 

To accomplish this, we make a few general observa- 
tions, which provide insight into the magnetic geometry 
and emission pattern. Firstly, we note that the shape of the 
main pulse changes significantly from an asymmetric single- 
peaked profile at 0.61 GHz to an asymmetric double-peaked 
profile at 1.414 GHz. This suggests that the magnetic geom- 
etry and possibly the beam pattern change with emission 
altitude. Both pulse profiles, however, are consistent with 
hollow cone emission. Secondly, in the main pulse, the lin- 
ear polarization follows the total intensity closely, suggesting 
that L oc cos S is a reasonable approximation. This is also 
true to a lesser degree in the interpulse. 

The stellar surface of this object lies at 0.13r-LC, placing 
a lower limit on the emission altitude. If we assume arbitrar- 
ily that the data at 0.61 GHz are already in the canonical 
basis, we find that, at (a, i) = (20°, 80°), the models for both 
pure and current-modified dipoles exhibit tilted balloons in 
the I-U and Q-U planes, similar to the data at r ~ 0.4rLC- 
Unfortunately, without additional information on the abso- 
lute orientation of f2p, we are limited to this assumption. 

Figures [9] and [To] show the pulse profile, PA swing, and 
Stokes phase portraits of (a) the main pulse and (b) the 
interpulse at 0.610 GHz and 1.414 GHz respectively. For the 
0.610 GHz case, we assume /3 = 0, whereas for the 1.414 GHz 
case, we align the narrow balloon shape in the Q-U plane of 
Figure [TOl a') with that of Figure [9{ a) by assuming j3 = 45°. 
In the top left panel of each subfigure, we plot ///max (solid 
curve), L//niax (dashed curve) and the PA swing (dotted 
curve) wherever L ^ O.ll/max. Stepping clockwise, the next 
three panels show I-Q, Q-U and I-U . 

We now examine the magnetic geometry, beam pattern, 
orientation, and emission altitude in more depth in Sections 
I4.1H4.4I We find that the model with a hollow cone and 
L (X cos 6 must be generalized by letting / and L vary with 
(f) in order to fit the data in detail. 



4.1 Magnetic geometry 

At 0.61 GHz, the Stokes phase portraits for the main pulse 
are all narrow balloons. The major axes of the balloons tilt 
in different directions: we find dQ/dl < in the I-Q plane, 
dU/dl < in the I-U plane, and dU/dQ > in the Q-U 
plane. Assuming /3 = 0, the tilt of the Q-U phase portrait 
discounts a pure dipole magnetosphere at a low emission 
altitude. At higher altitudes, where aberration is important, 
e.g. at r = O.lrLC, the pure and current-modified dipoles 
produce phase portraits that are similar enough at some 
orientations to warrant considering both cases. 

We now seek a match from the look-up tables for a 
hollow cone with L — I cos 9 (Figures I26H28|I , keeping in 
mind that we are interested in orientations which provide an 
interpulse (a or i > 70°). We find that, at (a, i) = (20°, 80°), 
the phase portraits for I-Q and I-U match approximately 
the balloons in the data, although the sign of U is reversed 
(this choice of orientation is justified in Section I4.2|l . 

4.2 Orientation (a,z) 

Finding the exact orientation is an iterative process, requir- 
ing the beam and polarization patterns to be adjusted at 
each step. Initially, we seek a match to the data at 0.61 GHz 
(Figure [9}. The interpulse is extremely useful in narrowing 
the range of possible orientations to 70° 90° (as- 

suming beams of intrinsically equal luminosity). As the in- 
terpulse is weaker than the main pulse, we know that a and 
i are less than 90° . From Figures I26H28I there are two ori- 
entations with similar balloons in all three phase portraits, 
namely {a,i) = (20°, 80°). The phase portraits for a pure 
dipole at {a,i) = (20°, 80°) are also similar. 

Before zooming in to refine the grid around (a, i) = 
(20°, 80°), we experiment with various emission altitudes 
while tailoring the pattern to fit the data. In Section 14.31 
we construct beam and linear polarization patterns at r = 
0.4rLC for the data at 0.61 GHz. 

4.3 Beam and polarization patterns 

Both the main pulse and interpulse at 0.61 GHz are single- 
peaked and skewed to the left. The interpulse peaks at 
~ 0.6/max. To capture this behaviour, we model the emission 
region as two hollow cones whose brightness varies longitu- 
dinally, i.e. the cones are shaped like horseshoes in cross- 
section. The best-fit beam pattern is given empirically by 

I{e,(f>) = (27rcri)-^/^ [0.8 -H|sin(0- 1.05)1] 

xexp [-0.5(9 - pif/af] (10) 

-f0.09(27rcr2)"^''^ [0.4 -f |sin((?!) - 1.65)|] 

xexp[-0.5(e-7r + p2)V'^2] , (H) 

where cti = 3° and a2 = 3.5° are the widths of the main 
pulse and interpulse respectively, and pi = 23° and p2 = 35° 
are the corresponding opening angles. The modelled pulses 
are « 3 times wider than the data as a result of the compu- 
tationally limited resolution of our numerical grid. 

Asymmetric emission regions are consistent with the 
patchy b eam model introduced to ex plain asymmetric pulse 
profiles ijLvne fc Manchester! Il988l ) and with theo retical 
models of pulsar magnetospheres like the slot gap (jAronsI 
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Figure 9. Polarimetry of (a) the main pulse and (b) the interpulse of PSR J1939+2134 at 0.61 GHz dStairs et al.lll999h . Each subfieure 
shows (clockwise from top left panel): (i) ///max (lower subpanel, solid curve) and /.//max (lower subpanel, dashed curve) profiles and 
PA swing (upper subpanel, dotted curve, in rad) versus pulse phase (in degrees); (ii) I-Q phase portrait; (iii) Q-U phase portrait; (iv) 
I-U phase portrait. Data are presented courtesy of the EPN. 
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Figure 10. As for Figure|9]but at 1.414 GHz IStairs et"allll999l) . Data are presented courtesy of the EPN 



1 19831 ). There are several successful precedents for pulse mod- 
els withJiorsesh££_beams^ the em pirical models proposed 
bv lKarastergiou fc JohnstonI (|2007l '). 

The linear polarization profiles of the main pulse and 
interpulse at 0.61 GHz look surprisingly different, naively 
suggesting a north-south asymmetry. In the main pulse, L 
follows the pulse profile closely, lagging the pulse centroid in 
phase by « 4.5°, peaking at « 0.87max. In the interpulse, L 
is extremely low, peaking at « 0.1/max, and appears to be 
triple-peaked. Despite the apparent difference in the profiles, 
we are able to reproduce them surprisingly well using the 
same model, given by 

L(6l,(^) = |cos6lsin(<^-(-0.92)| (12) 
without invoking a north-south asymmetry. As expected. 



however, p2p reproduces the L profile of the main pulse 
more accurately than that of the interpulse. We emphasize 
that (|ll|l and (|12|) are certainly not unique and do not fit 
the data exactly, but they are adequate for the empirical 
task at hand. 

Adopting (lllf) and (|12|) . we generate zoomed-in look- 
up tables for both pure and current-modified dipoles, in the 
range 14° ^ a sC 24°, 76° i ^ 84°, with a resolution of 2°. 
We find the closest match is for a current-modified dipole 
at (a, i) = (22°, 80°), with a 'by eye' uncertainty of ±2° for 
a and ±1° for i. This margin would widen if I{8,(j)) and 
L{6, <j>) were adjusted for each orientation. 

In Figure 1111 we plot the pulse profile, PA swing and 
Stokes phase portraits of the model at r = 0.4rLc and 
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(a,j) — (22°, 80°). The slight jaggedness of the pulse pro- 
files is a product of the finite grid resolution. The Stokes 



phase portraits of the main pulse [Figure If (a) match the 



data in Figure 9(a) reasonably well. In the data, the I-Q 
balloon ranges from —0.5 ^ Q Q, whereas in the model 
it is thinner and ranges from —0.3 Q 0. The I-U bal- 
loon in the data ranges from —0.6 li^U <Q, whereas in the 
model it ranges from —0.65 "^^U <Q. The PA swing in the 
data is nearly flat, with a slight negative gradient, whereas 
the model shows a slight positive gradient. For the inter- 
pulse, there is poorer agreement in L. The tilted balloon 



in I-Q from the data [Figure 9(b) is reproduced in Figure 
ll(b)[ including the kink visible at (/, Q) ^ (0.2,-0.03). 



In the I-U plane, the data feature a tilted balloon, with a 
prominent kink at (/, U) ~ (0.2, —0.02). Our model shows a 
broad hockey stick instead. At a stretch, one may perhaps 
argue that the balloon with the kink resembles the hockey 
stick qualitatively, but we do not press the point. We note 
that if the upper half of the hockey stick is reflected about 
(7 = 0, it would match more closely. In the Q-U plane, the 
tilted oval seen in the data is reproduced in the model, but 
with -0.07 <U < 0.05 in the data, and -0.08 < (7 < in 
the model. The modelled PA swing is flat, with a negative 
gradient, and lacks the phase- wrapping seen in the data. 

4.4 Emission altitude 

According to the standard radius-to-frequency mapping, the 
observation fr equency scales with emission radius as r~^^'^ 
(|Ruderman fc Sutherland 1975; Cordcs 1978). If the data at 
0.61 GHz correspond to r = 0.4rLC, then 1.414 GHz corre- 
sponds to r = 0.22rLC. Figures [T2I413I show the pulse pro- 
files, PA swing and Stokes phase portraits predicted theo- 
retically for both pulses, for emission altitudes ranging from 
r = 0.22rLC to 0.35rLC. The relative heights of the pulses 
change with emission altitude. We label them Pulse 1 (Fig- 
ure [T^]), corresponding to the main pulse in the data, and 
Pulse 2 (Figure [13)) . corresponding to the interpulse in the 
data. 

The theoretical pulse profile and phase portraits at 
r = 0.22rLC (Figures I12H13I top row) display some inter- 
esting features. First, the main pulse and interpulse have 
roughly the correct shapes, but swap positions in phase, i.e. 
the hollow cone which emits the main pulse at r = 0.4rLC 
also emits the interpulse at r = 0.22rLC, and vice versa. 
Upon inspection, it is likely that the same is true in the 
data. The triple-peaked linear polarization profile seen in 
Figure 9(b) is also present in Figure 10(a) although the 



first component in L is much weaker than the second and 
third at 1.414 GHz. Additionally, the kinks seen in the I-Q 
and I-U planes of Figure 9(b) are seen in Figure 10(a) In 
the data, the interpulse peaks at « 0.5/max, compared to 
K, 0.97inax in the model at r = 0.22rLC. 

Second, the linear polarization profile and phase por- 
traits at r = 0.22rLC reproduce the main pulse reasonably 
well but do not match the interpulse. The linear polariza- 
tion of the main pulse (Pulse 1; Figure [T^ top row) is ~ 50% 
weaker than observed [Figure |10(a)| . Also, in the simulated 
profile, the third L peak in the main pulse is compara- 
ble in height to the second peak, while in the data it is 
weaker. In the simulated I-Q and I-U planes of the main 
pulse, there are reasonable matches to the balloons in the 



data. In the simulated I-Q plane, we see a balloon with 
a kink at {I,Q) ~ (0.6,-0.1). In the data, the kink ap- 
pears at {I,Q) ~ (0.4,-0.1). In the I-U plane, the kink 
seen in the data at {I,U) ~ (0.4,-0.1) is reproduced at 
(7, U) k: (0.6, -0.2) in the model. In the Q-U plane, the 
data trace out a thin balloon with dU/dQ > spanning 
-0.3 < Q < and -0.4 <U < 0.05. The simulated phase 
portrait shows a thin balloon with the same orientation, 
spanning -0.15 < Q < and -0.3 < U < Q. For the in- 
terpulse (Pulse 2; Figure 1131 top row) the simulated total 
intensity is twice the observed intensity, and the simulated 
linear polarization is « 6 times stronger than observed [Fig- 
The simulated and observed balloons in the I-Q 



10(b) 



and I-U planes are rotated by 90° clockwise with respect to 
the data, whereas the Q-U balloon is rotated by 180°. These 
discrepancies are also refiected in the PA swing. 

As the emission altitude increases from r = 0.22rLC to 
r = 0.35r-LC, the phase portraits of the main pulse (Pulse 
1; Figure [T^J change. In the I-Q plane, the kink in the bal- 
loon shifts towards (/, Q) ^ (0.3,-0.08). At r = 0.35rLC, 
the balloon in the I-U plane starts to resemble the hockey 
stick seen in Figure [11(b) | In Q-U, the thin balloon rotates 
clockwise. In the interpulse (Pulse 2; Figure [13}, the I-Q 
balloon narrows and lengthens in Q, and the Q-U balloon 
narrows. 



From Figures rr2H13l we draw the following conclusions. 

(i) Although the simple model given by and (|12|) mod- 
els the 0.61 GHz data reasonably successfully, it fails for the 
data at 1.414 GHz displayed in Figure [10] However, the ob- 
served pulse profiles and Stokes phase portraits suggest that 
the emission region of the main pulse at 0.61 GHz corre- 
sponds to that of the interpulse at 1.414 GHz, and vice versa. 
Additionally, the emission pattern may change with radius. 

(ii) The discrepancies between the data and the phase por- 
traits at r = 0.22rLc (the altitude predicted by the radius- 
to-frequency mapping) indicate that the toroidal field may 
not increase monotonically with r. The phase portraits for 
the interpulse between r = 0.22rLC and r — 0.35rLC are all 
a poor match to the data, (iii) It is possible that the data 
should be referred to a different value of /3 at 1.414 GHz 
than the one we assume, which would rotate the Q-U phase 
portrait, and change the shapes of the I-Q and I-U patterns. 



5 A DETAILED MULTI-PEAKED EXAMPLE: 
PSR J0437-4715 

We now repeat the procedure in Section |¥] and CMIO for 
PSR J0437-4715. This object has P = 5.758 ms and P = 
5.729 X 10"^° s s"^ ([Bell et al.lll997l '). It was chosen because 
it exhibits five distinct peaks in its pulse profile, clearly vis- 
ible at 1.44 GHz, and a highly structured PA swing. There 
is no interpulse observed in this object. Unlike the other 
objects considered in this paper and CMIO, we find that 
PSR J0437— 4715 cannot be modelled by either a pure or a 
current-modified dipole field, even if a multiple-peaked beam 
pattern is constructed emiprically to fit the I{t) data ex- 
actly. Indeed, the Stokes phase portraits point persuasively 
to the existence of a strong quadrupole and higher-order 
multipoles at the radio emission altitude. In this respect, 
PSR J0437— 4715 is an excellent candidate for more detailed 
Stokes tomography studies in the future. In this section, we 
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Figure 11. Theoretical polarization model of (a) the main pulse and (b) the interpulse of PSR J1939+2134 for a current-modified dipole 
emitting at r = 0.4rLCi with orientation (a, i) = (22°, 80°), beam pattern given by lllll l. and linear polarization given by I I12I I. Each 
subfigure shows (clockwise from top left panel): (i) ///max (lower subpanel, solid curve) and L//max (lower subpanel, dashed curve) 
profiles and PA swing (upper subpanel, dotted curve, in rad) versus pulse phase (in degrees); (ii) I-Q phase portrait; (iii) Q-U phase 
portrait; (iv) I-U phase portrait. 



restrict ourselves to presenting the argument that the pure 
and current-modified dipoles categorically fail to match the 
data for the polarization models that work well for the other 
objects studied in this paper and CMIO. 

Data for PSR J0437-4715, at 1.44 GHz and 4.6 GHz, 
are obtained from the EPN online archive. The data were 
originally published bv lManchester fc JohnstonI (|l995l ). Fig- 
ure [14] presents the pulse profile, phase portraits, and PA 
swing (in rad) at 1.4 GHz. The top left panel shows ///max 
(solid curve), L//max (dashed curve) and the PA swing (dot- 
ted curve, at longitudes where L ^ O.lLmax). Stepping 
clockwise, the next three panels show I-Q, I-U and Q-U.^ 

As with PSR 1939-1-2134, we assume that the data ob- 
tained from the EPN are presented in the canonical polar- 
ization basis at one reference frequency, chosen here to be 
1.44 GHz. The stellar surface is at r = 0.036r-LC, providing 
a lower limit on the emission altitude. 

The pulse profile has five components, labelled A-E 
in Figure 1141 two on either side of the largest peak (C). 
Each component in the pulse profile corresponds to a dis- 
tinct sub-pattern in the Stokes phase portraits. In the I-Q 
plane, peaks A and B correspond to the two small loops at 
Q > 0, peak C is the large figure-eight, and peaks D and E 
correspond to the loop at Q < 0. In the I-U plane, peaks 
A and B correspond to the kinks at (7 < 0, peak C is the 
large balloon, and peaks D and E correspond to the kink 
at [/ > 0. The U-Q plane is complicated, forming a rough 
X-shape, with one diagonal having dU /dQ > (we call this 
diagonal 1), and the other having dU/dQ < (diagonal 2). 
Diagonal 1 corresponds to peak C. In diagonal 2, peaks A 
and B occupy the Q > region, whereas peaks D and E 
occupy the Q < region. 
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Figure 14. Polarimetry of PSR J0437-4715 at 1.44 GHz 
^Manchester &: JohnstonI Il995h . Clockwise from top left panel: 
(a) ///max (lower subpanel, solid curve) and ////max (lower sub- 
panel, dashed curve) profiles, and PA swing (upper subpanel, dot- 
ted curve, in rad) all plotted against pulse phase (in degrees); (b) 
I-Q phase portrait; (c) Q-U phase portrait; (d) I-U phase por- 
trait. Data are presented courtesy of the EPN. 



The linear polarization within peak C is double-peaked. 
This kind of structure is common and is modelled adequately 
by a filled core beam with L = I sin 6, as demonstrated for 
several objects in CMIO. The phase separation of the peaks 
suggests that A, B, D, and E originate from two hollow cones 
centred on the same axis (peak C). Peak pairs B/C and C/D 
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Figure 12. Theoretical polarization model of Pulse 1 of PSR J1939+2134 as a function of emission altitude, for a current-modified 
dipole with (a,i) = (22°, 80°). In landscape mode, the plots for each emission altitude occupy rows, increasing from r = 0.22rLc (top 
row) to r = 0.35rLC (bottom row). From left to right, the columns show (1) ///max (solid curve) and L//max (dashed curve) versus 
pulse phase I (in units of degrees), (2) I-Q phase portrait, (3) I-U phase portrait, (4) Q-U phase portrait, and (5) the PA swing (in rad; 
data points with L ^ O.lLmax plotted only). 
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Figure 13. As for Figure [12] but for Pulse 2 of PSR J1939+2134. 



are separated by ~ 0.7 rad and ~ 0.9 rad respectively, while 
A/C and C/E are both separated by ~ 1.2 rad. We confirm a 
posteriori that a and i lie in ranges where interpulse emission 
does not contribute significantly. 



5.1 Magnetic geometry 

In order to determine the magnetic geometry, we search the 
look-up tables for a good match involving a filled core and 
a hollow cone beam. We assume L = I sm9 because of the 
double-peaked L profile in peak C. The phase portraits for 
the filled core should match the large patterns corresponding 
to peak C, while the phase portraits for the hollow cone 
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should match the smaller patterns. We do not expect perfect 
matches due to the complex beam and polarization patterns. 
At a minimum, however, we seek an approximate match 
for the rough figure-eight that forms diagonal 1 in the Q-U 
plane, the figure-eight in I-Q, and the balloon in I-U. 

The pure dipole look-up tables at r = O.frLC do not 
feature figure-eight shapes in the I-Q plane for any orienta- 
tion. For a < i, there are balloons in I-U, and heart shapes 
in Q-U. As in many other opbjects, a pure dipole is ruled 
out (CMIO). 

For the current-modified dipole (Figures l22H24p . the 
closest match is at {a,i) = (30°, 30°) (see Section [5.21 for a 
detailed justification). For a filled core, there are asymmetric 
figure-eights in the Q-U and I-U planes fFigures l221424p . and 
a broad oval in I-Q (Figure [22]). The hollow cone phase por- 
traits at this orientation feature asymmetric mosquitoes in 
I-Q and I-U (Figure [30[431|) and an asymmetric heart shape 
in Q-U (Figure [32|) . At this stage, we cannot confidently dis- 
count either magnetic configuration, as the Stokes portraits 
change when the filled core and hollow cones are combined. 
This issue is examined thoroughly in Sections 15. 3H5. 41 



5.2 Orientation 

In this section, we justify {a,i) = (30°, 30°) as the best 
matching orientation. As PSR J0437— 4715 does not have 
an observed interpulse, we rule out all orientations with 
Q ^ 80° , i > 80° . We also rule out orientations with i > a 
because the associated phase portraits look nothing like the 
data. For example, for a filled core, Figures I22H24I contain 
balloons in I-Q, narrow, tilted balloons and straight lines in 
I-U, and ovals in U-Q. None of these patterns appear in the 
data in Figure [HI The Q-U discrepancy is especially signif- 
icant as the shape of the Q-U portrait is basis-independent. 

The best matching orientation predicted by the R VM 
is (a,i) = (145°, 140°) jManchester fc Johnstonll 19951 ). al- 
though the authors note that the PA swing deviates largely 
from the model. For this reason, no formal uncertainties are 
assigned to the fitted parameters, which were chosen to be 
a reasonable representation of the data. For a dipole field, 
the Stokes phase portraits are symmetric about (q, i) = 
(90°, 90°), i.e. the phase portraits for (a,i) = (145°, 140°) 
and (35°, 40°) are identical. From CMIO, the phase portraits 
for a pure dipole at (q, i) — (40°, 40°) and r ^ tlc, with 
a hollow cone and L — sin 9, feature a narrow, tilted bal- 
loon in I-Q, a mosquito in I-U, and a heart in Q-U (see 
Figures 22-24 in CMIO). The I-U and Q-U shapes are sym- 
metric about U = 0. Interestingly, although we have chosen 
our best match independently of the RVM results, the two 
orientations are close. 

We reiterate that the complex multiple-peaked beam 
and polarization patterns complicate the matching process. 
Some orientations must be tested with beam patterns tai- 
lored to fit the data, as described in Section r5.3l before being 
ruled out. For (q, i) = (30°, 30°), the appropriate beam pat- 
tern is a filled core surrounded by two hollow cones. The 
resulting phase portraits show distorted, tilted balloons in 
both I-Q and I-U. These shapes resemble roughly the data 
in Figure 1141 although there are large discrepancies too, 
chiefiy that the figure-eight in I-Q is missing, and that I-U 
is not symmetric about U — 0. In Sections 15.31 and 15.41 we 



construct detailed beam and linear polarization models in 
an attempt to improve the fits. 

5.3 Beam pattern 

In fitting the complex puls e profile o f PSR J0437-4715 
at 1.44 GHz, Gangadhara fc Thomad (j2008l ) identified 11 
Gaussian components. They proposed that the beam pat- 
tern comprises five nested cones at different altitudes cen- 
tred on the filled core, and that the altitudes range from 
0.07rLC-0.3rLC. 

We consider a simpler model and focus on one fixed 
altitude. We model the pulse profile empirically with a filled 
core, Ii{6,(j)) (peak C), surrounded by two hollow cones, 
l2{9, (t>) (peaks B and D) and hid, 0) (peaks A and E). The 
intensity maps take the empirical forms 

h{e,^) = (2^a?)-i/2 {exp [-eV(2a?)] 

+exp[-(e-^)V(2a?)]}, (13) 
l2{e,4>) = /32(<^)(2^a2')-^/2{cxp[-(e-p2)V(2(7i)] 

+e^p[-{e-7v + p2f/i2al)]}, (14) 
hie, CP) = /33ic^){27valr'^^{e^p[-ie-p;f/{2al)] 

+exp[-{e-Tv + p:if/i2al)]}, (15) 

where ai — 6.5°, (T2 = 2.5° and as — 2° are the beam widths 
of the core and cones, p2 = 18° and pa — 26° are the opening 
angles of the two cones, and l32{4>) and P^iij)) are functions 
describing the longitudinal structure of the two cones, given 
empirically by 

/32((^) = 0.06|cos(0.3(^)|, (16) 
/33(<^) = 0.06|cos(0.75(^)|. (17) 

As in Section |4l the cones are shaped like horseshoes. Given 
(|13|) - (|15p . we also find that the linear polarization pattern 
is fitted empirically by 

Lie,(j}) = (361)"^ sin(6l- 0.01) +0.07|cos(0.75<|i)| 

X (2^a|)-^/^exp [-(S - P2)V(2<x|)] • (18) 

We emphasize that equations (|13|) - (|18p are not unique 
fits, nor do they produce perfect agreement with the data. 
In particular, the data show that the B/C peaks are closer 
to each other than C/D, yet we are unable to reproduce 
this with a reasonably simple model. Our modelled peaks 
are equidistant. The models are sensitive to the pulsar's 
orientation. Every time we vary a or i around (30°, 30°), 
we must adjust the coefficients in (|13[) " (|18p . We find that 
the closest match to the data, although poor, is achieved at 
(a,j) — (32°, 26°). We note as well that the Stokes phase 
portraits match marginally better if we rotate the polariza- 
tion basis by /? = 90°. The difficulty in achieving a good 
match may well be telling us that the underlying magnetic 
geometry is not a current-modified dipole. 

5.4 Decomposed phase portraits 

We now demonstrate how the Stokes phase portraits change 
as we add h and I3 to the filled core Ji. Figure [151 shows the 
pulse profile, PA swing, and phase portraits for h, Ii 12, 
and h + I2 + h respectively at an emission altitude of r = 
O.lrLC and /3 = 90°. 
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For just the filled core (Figure [151 top row in landscape 
orientation), corresponding to peak C in the data, there is 
a thin, tilted balloon in I-Q, a balloon in I-U , and a tilted 
oval in Q-U . Aside from the obvious dissimilarity with the 
figure-eight seen in the data in I-Q, the oval in Q-U does not 
resemble diagonal 1 in the data. The balloon in I-U is an 
approximate match to the data, although it is not symmetric 
about U = Q. 

For the filled core and one hollow cone (Figure fTSl mid- 
dle row), i.e. peaks B, C and D in the data, the conal com- 
ponents introduce a kink at (/, Q) ~ (0.2,-0.03). In the 
I-U plane, kinks are predicted to occur at / < 0.2. The data 
also contain kinks in this region. In the Q-U plane, there is 
another kink near (Q, U) ~ (0, —0.05), though the model is 
still a poor match to the data. 

The addition of the second hollow cone in Figure [T5] 
(third row) completes the beam pattern. Still, the overall 
shapes of the phase portraits do not match the data. The 
second cone appears as an extra kink in I-Q and I-U in 
the region / < 0.15. In Q-U, a secondary oval forms. At a 
stretch, it might be said that this secondary oval corresponds 
to one diagonal of an X-shape while the other large oval 
corresponds to another, but other interpretations are equally 
possible. Adjusting the emission altitude does not improve 
the fit. 

For completeness, in Figure 1161 we present the pulse 
profile, PA swing, and phase portraits for a filled core with 
two hollow cones for a pure dipole field at the same ori- 
entation and altitude, and /3 = 0. The phase portraits are 
also a poor match to the data. The I-Q plane features a 
tilted, asymmetric balloon, whereas the I-U plane features 
an asymmetric figure-eight. In the Q-U plane, there is a dis- 
torted oval surrounded by a tilted heart shape. Again, the 
basis-independent Q-U shape does not resemble the data at 
all. The kinks seen in the phase portraits of the current- 
modified dipole also appear in the pure dipole. 

Finally, in Figure [17] we present phase portraits for a 
simplified version of the current-modified field where the 
toroidal field is given by 



-B<4 ~ 



(19) 



In this stripped-down expression, depends on 6 and a 
only through the poloidal field, and scales simply as r/rLC- 
The phase portraits in Figure [17] are also presented for a 
filled core with two hollow cones at the same orientation 
and altitude, and /3 = 0. Again, they are a poor match to 
the data, although the I-Q plane now features a tilted figure- 
eight similar to the data. The I-U plane features an asym- 
metric figure-eight, and the Q-U plane features two inter- 
locking ovals. As in the previous cases, kinks corresponding 
to the various pulse peaks punctuate the phase portraits. 

We note that, for the three magnetic geometries con- 
sidered, the theoretical PA swings are smooth and relatively 
flat and do not contain any of the kinks seen in the data. 
Attempts to rotate the Q-U portraits of the pure and simpli- 
fied current-modified dipole to yield a better fit, i.e. trialling 
several values of /3, are also unsuccessful. For example, for 
p = 150° with the simplified current-modified dipole, diag- 
onal 1 in the Q-U plane aligns with the large, interlocking 
ovals in the model while diagonal 2 aligns with the smaller, 
third oval. However the I-Q plane now features a large bal- 
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Figure 16. Theoretical polarization model of PSR J0437— 4715 
for a pure dipole emitting at r = O.lrLc with {ct,i) = (32°, 26°), 
beam pattern given by I I13I I. I I14I I and JTSj, and linear polarization 
given by I I18II . Clockwise from top left panel: (a) ///max (lower 
subpanel, solid curve) and L (lower subpanel, dashed curve) pro- 
files, and PA swing (upper subpanel, dotted curve, in rad) all 
plotted against pulse longitude (in units of degrees) (b) I-Q phase 
portrait; (c) Q-U phase portrait; (d) I-U phase portrait. 



loon symmetric about Q = 0, and the I-U plane features a 
broad, tilted balloon, neither of which matches the data. 

We conclude that we are unable to fit the Stokes phase 
portraits for PSR J0437— 4715 satisfactorily using a pure 
or current-modified dipole and a wide range of trial-and- 
error models for the beam and linear polarization pat- 
terns. There are several possible reasons for this, (i) The 
beam and polarization patterns might be very different, 
e.g. two highly asymmetric, nested cones, (ii) If the ob- 
served emission does in deed originate from multiple alti- 
tudes (jCangadhara fc Th omas 2008), we would be unable 
to reproduce the Stokes phase portraits even if our guesses 
for a,i,I{t) and L{t) are correct, (iii) The magnetic field in 
the emission region is neither a pure nor a current-modified 
dipole (very likely). 



6 CONCLUSION 

In this paper, we generalize the Stokes tomography tech- 
nique introduced by CMIO by adding interpulse emission. In 
Section|3] we present the Stokes phase portraits of 15 MSPs 
from the EPN online archive. By comparing the data to the 
generalized look-up tables for a current-modified dipole in 
the Appendix, we are able to infer approximately the ge- 
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Figure 15. Bottom-up, component-wise assembly of theoretical polarization models for PSR J0437— 4715 for a current-modified dipole 
emitting at r = O.lrLc with {a,i) = (32°, 26°) and /3 = 90°. In landscape orientation, from top to bottom, we plot the following beam 
patterns: (top row) peak C, filled core, equation (middle row) peaks B-D, filled core plus hollow cone, equations and II14I I: and 

(bottom row) peaks A-E, filled core plus two hollow cones, equations II13I I. II14I I and II15I I. Linear polarization in all three rows is given 
by JTHJ. From left to right, in landscape orientation, the columns contain (1) ///max (solid curve) axid L/Xmax (dashed curve) profiles. 
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Figure 17. Theoretical polarization model of PSR J0437-4715 
for a simplified current-modified dipole given by I I19I I emitting at 
r = O.lrLC with (a, i) = (32°, 26°), beam pattern given by I I13I I. 
II14I I and I I15II . and linear polarization given by I I18I I. Clockwise 
from top left panel: (a) ///max (lower subpanel, solid curve) and 
L (lower subpanel, dashed curve) profiles, and PA swing (upper 
subpanel, dotted curve, in rad) all plotted against pulse longi- 
tude (in units of degrees) (b) I-Q phase portrait; (c) Q-U phase 
portrait; (d) I-U phase portrait. 



netic mountain (|Pavne fc Melato3l2004l : IVigelius fc Melatoj 
[2008). 

We emphasize the utility of the Stokes phase portraits 
as a supplementary diagnostic tool for MSPs. The PA swing 
on its own is especially ambiguous when dealing with non- 
dipolar fields. Future work will focus on the role played by 
circular polarization in Stokes tomography, the longitudinal 
structure of vacuum and force-free magnetospheres, and the 
polarization signatures of magnetic mountains. These topics 
will form the subject of companion papers. 
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ometric orientations for five of the MSPs. This is an im- 
provement on the PA swing and rotating vector model, 
which yield orientations for only two of the objects — 
orientations which, it transpires, are inconsistent with the 
observed Stokes phase portraits. In Section |4l we model 
PSR J19394-2134 in detail, obtaining a match for the data 
at 0.61 GHz with a current-modified dipole for (a, i) — 
(22± 2°, 80± 1°) and r = 0.4rLC- However, we are unable to 
reproduce the data at 1.414 GHz for the same orientation at 
altitudes in the range 0.22r-LC =^ ^ 0.35rLc- In Section [SI 
we repeat the process with PSR J0437-4715. At 1.44 GHz, 
even the closest-matching orientation, {a,i) = (32°, 26°) 
with r = O.lrLC, does not reproduce the data satisfactorily. 

The results from Sections [4] and [5] indicate that, while 
pure or current-modified dipoles are effective models for 
non-recycled pulsars (CMIO), MSPs are likely to have more 
complicated magnetic geometries. This is not surprising, 
as the accretion p rocess can significantly distort a pul- 
sar's mag netic field (|Lai et al ] |l999l : iPavne fc Melatosll2004l : 
iLamb et a l. 2009). Alternative magnetic configur ations in- 
clude a quadrupole or localized surface a nomaly (iLai et al.l 
19991: iLong et al."2008'). a force-free field (|Spitk ovskv"2006|; 



Bai fc Sp itkovskv 2009), a vacuum- like field ((Melatos 1993), 



or a field distorted by the formation of a polar mag- 
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APPENDIX: ATLAS OF LOOK-UP TABLES OF 
STOKES PHASE PORTRAITS 

We present look-up tables for the beam patterns and lin- 
ear polarization models used in CMIO, updated to include 
interpulse and relativistic aberration effects. All figures are 
for r = O.lrLC- Stokes phase portraits and PA swings are 
shown for a current-modified dipole with 

(i) a filled core beam with L — I cos 9 (Figures If 8H2ip . 

(ii) a filled core beam with L = I sin 9 (Figures I22H25|I , 

(iii) a hollow cone with L — I cos 9 (Figures I26H29|I , and 

(iv) a hollow cone with L = I sin 9 (Figures I30H33P . 
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Figure 18. Current-modified dipole. Look-up table o: 
polarization L = I cos 8, where is the emission point 
of increasing 10° ^ i 90° (left-right) and 10° ^ a 
normalised by its peak value. Q is plotted on the ver 
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Stokes phase portraits in the I-Q plane for filled core beams with degree of linear 
colatitude, and r = O.lrLC- The panels are organised in landscape mode, in order 
f ^ 90° (top-bottom) in intervals of 10°. I is plotted on the horizontal axis and 
ical axis. 
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Figure 19. Current-modified dipole. Layout as for Figure [181 but for I-U (I on the horizontal axis). 
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Figure 20. Current-modified dipole. Layout as for Figure (TS] but for Q-U {Q on tfie horizontal axis). 
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Figure 21. Current-modified dipole. Layout as for Figure [Tsl but for position angle (on the vertical axis in landscape orientation, 
units of radians) versus pulse longitude (on the horizontal axis, in units of 2-k radians). 
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Figure 22. Current-modified dipole. Look-up table of Stokes phase portraits in the I-Q plane for filled core beams with degree of linear 
polarization L = 1 sinS, where 9 is the emission point colatitude, and r = O.lrLC- The panels are organised in landscape mode, in order 
of increasing 10° ^ i ^5 90° (left-right) and 10° ^ a ^ 90° (top— bottom) in intervals of 10°. I is plotted on the horizontal axis and 
normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 23. Current-modified dipole. Layout as for Figure [22l but for 1-U (I on the horizontal axis). 
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Figure 24. Current-modified dipole. Layout as for Figure [22] but for Q-U {Q on tfie horizontal axis). 
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Figure 25. Current-modified dipole. Layout as for Figure [22l but for position angle (on the vertical axis in landscape orientation, 
units of radians) versus pulse longitude (on the horizontal axis, in units of 2-k radians). 




Figure 26. Current-modified dipole. Loolc-up table of Stolces piiase portraits in tiie 1-Q jjlane for liollow cones with opening angle 25° 
and degree of linear polarization L = /cos 9, where 9 is the emission point colatitude, and r = O.lrLC- The panels are organised in 
landscape mode, in order of increasing 10° ^ « ^ 90° (left-right) and 10° ^ a ^ 90° (top— bottom) in intervals of 10°. I is plotted on the 
horizontal axis and normalised by its peak value. Q is plotted on the vertical axis. 



C. T. Y. Chung et al. 



diti^ 



\ 



!-i r= p - 



i ^. 



ill. 



...-■1. 



I ■ 

J: 



./I 



Figure 27. Current-modified dipole. Layout as for Figure [26l but for I-U (I on the horizontal axis). 
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Figure 28. Current-modified dipole. Layout as for Figure [26l but for Q-U {Q on tfie horizontal axis). 
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Figure 29. Current-modified dipole. Layout as for Figure [26l but for position angle (on the vertical axis in landscape orientation, 
units of radians) versus pulse longitude (on the horizontal axis, in units of 2-k radians). 
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Figure 30. Current-modified dipole. Loolc-up table of Stolces piiase portraits in tlie I-Q plane for hollow eones with opening angle 25° 
and degree of linear polarization L = Isind, where is the emission point colatitude, and r = CItlc. The panels are organised in 
landscape mode, in order of increasing 10° C i ^ 90° (left-right) and 10° ^ o ^ 90° (top— bottom) in intervals of 10°. / is plotted on the 
horizontal axis and normalised by its peak value. Q is plotted on the vertical axis. 
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Figure 31. Current-modified dipole. Layout as for Figure [30l but for I-U (I on the horizontal axis). 
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Figure 32. Current-modified dipole. Layout as for Figure l30l but for Q-U {Q on tfie horizontal axis). 
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Figure 33. Current-modified dipole. Layout as for Figure [30l but for position angle (on the vertical axis in landscape orientation, 
units of radians) versus pulse longitude (on the horizontal axis, in units of 2-k radians). 
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